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Abstract 

A study of the leeside flow characteristics of the 
Shuttle Orbiter is presented for a reentry flight con- 
dition. The flow is computed using a point-implicitT 
finite-volume scheme known as the Langley Aerother- 
modynamic Upwind Relaxation Algorithm (LAURA). 
LAURA is a second-order accurater laminar Navier- 
Stokes solverTincorporating finite-rate chemistry with a 
radiative equilibrium wall temperature distribution and 
finite-rate wall catalysis. The resulting computational 
solution is analyzed in terms of salient flow features and 
the surface quantities are compared with flight data. 
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Acronyms 

AFRSI Advanced flexible reusable surface insulation 
DFI Developmental flight instrumentation 
LAURA Langley aerothermodynamic upwind 
relaxation algorithm 
OMS Orbital maneuvering system 

PNS Parabolized Navier-Stokes 

RCC Reinforced carbon-carbon 

RCG Reaction cured gas 

RCS Reaction control system 

SILTS Shuttle infrared leeside temperature sensing 
STEIN Supersonic/hypersonic three-dimensional 
external inviscid flow 
STS Space transportation system 

UTC Universal coordinated time 
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Symbols 

B Semi-span of the wingFm 

C Local chord of the wingFm 

h Altituderkm 

L Reference lengthrm 

Moo Mach number 

Poo Freestream pressureFPa 

q Heat fluxrW/cm 2 

qo Reference heat fluxrW /cm 2 

Q Freestream dynamic pressureFPa 

T TemperaturerK 

To Reference temperatureFK 

Voo Freestream velocityTm/s 

X Distance along the longitudinal axisFm 

Y Distance along the pitch axisFm 

a Angle of attackrdeg 

[3 Yaw anglerdeg 

e Emissivity 

T jv Atomic nitrogen recombination coefficient 
T o Atomic oxygen recombination coefficient 
Poo Freestream densityTkg/m 3 
a Stefan-Boltzmann constantrW/cm 2 K 4 
$ Meridional anglerdeg 

Introduction 

During design of lifting-body configurations such as 
the HL-20rHermesrHope HTBuranretcTthe harsh en- 
vironment during reentry dictates a significant challenge 
of coping with the associated high heating rates on the 
windward surfaces. As a resultTthe windward aerother- 
modynamic environment has been extensively studied 
and is fairly well understood. Correspondinglyrthe ther- 
mal protection system (e.^T reusable tiles) can be de- 
signed rather efficiently with regard to both safety and 
weight requirements since the windward environment is 
reasonably well-defined. The leeside environmentThow- 
everrhas not been studied in great depth. Whereas the 
extremely high heating values are not commonplace on 
the leeward side of the vehiclerthe flow physics remains 
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largely unexplored. Due to this vague knowledge of the 
leeside flowfieldrthe design parameters have to be very 
conservative. This implies a significant weight penalty; 
and as Bill LearT of Learjet. fame would sayr “I’d kill 
grandma for another pound of weight [reduction].” 

The inability of ground test facilities to reproduce 
the high energy flows present during reentry conditions]” 
coupled with the prohibitive expense of flight testsldeads 
to the use of analytic, methods to describe the flow. 
Several previous efforts have been made to compute 
complex vehicle flowfields for reentry speeds for vari- 
ous conditions. 1-1 b Founded on some of the principals 
of Kutler et al. I 1 the STEIN code of Marconi et a/.P was 
among the first to treat a Shuttle-like configuration using 
inviscid analysis. Shortly thereafteiTSchiff and St.eger 3 
introduced the subsonic sublayer approximation for the 
Parabolized Navier-Stokes (PNS). This enabled stable 
space-marching of the PNS equations along a body. The 
PNS work was continued by Li 4,5 which consisted of a 
reentry flow around the a Shuttle-like vehicler incorpo- 
rating “real” gas effects. Venkatapthy b used the PNS 
formulation to compare with wind tunnel data for the 
Orbiter. This work was followed by Szema et al. I 10 
who computed a reentry flight condition for a Shuttle- 
like configuration. Balakarishnan' presented a zero an- 
gle of attackrhigh Mach number conditionrincorporat- 
ing “real” gas effects. Prabhu et al. I 8 computed a reen- 
try flight condition for the Shuttle Orbiter and com- 
pared with flight data. Meanwhilerthe more computer- 
memory intensive Euler equation solvers where also be- 
ing developed and applied to Shuttle-like vehicles at 
high speeds. 9-13 More recentlyFsufficient computational 
power has enabled the solution of the Full or Thin-layer 
Navier-Stokes equations over complex reentry vehicles to 
be computedrincluding finite-rate chemistry. 14-1 * 5 

A limited amount of flight data is available for 
the Shuttle Orbiter. For the first five flights (STS-1- 
5) of the Shuttle Orbiter ColumbiaFthe vehicle was in- 
strumented to measure surface pressures and tempera- 
tures. This instrumentation is know as the Developmen- 
tal Flight Instrumentation (DFI). The surface temper- 
atures and related heating rates are given by Hartung 
et a/. 1 ' -19 In additionFon subsequent flights of Shuttle 
Orbiter ColumbiaFthe Shuttle Infrared Leeside Temper- 
ature Sensing (SILTS) experiment 20 provides a map of 
the surface temperatures over portions of the leeside sur- 
face of the vehicle. 

This paper intends to make comparisons with avail- 
able flight data for the Shuttle OrbiteiT and to provide 
some insight of the intricate leeside flow structures dur- 
ing reentry. This flowfield analysis serves to assess the 
ability of the computational models within LAURA to 
predict the leeside environment through examining not 
only the leeside pressures and heating ratesr but also 



Figure 1: Shuttle Orbiter surface definition. 



Figure 2: Simplified Shuttle Orbiter surfaeb 
definition. 


through a detailed comparison of the temperature distri- 
butions predicted for the leeside of the wing as compared 
with flight data. 

Shuttle Geometry 

The actual Shuttle Orbiter Geometry is shown in 
Fig. 1. Since the complexity of the aft portion repre- 
sents a very arduous geometric modeling problemFit was 
simplified for the purposes of this investigation. The geo- 
metrical simplifications consist of the omission of the tail 
surfacesrbody flapsFand continuing the wing’s hinge-line 
edge thickness as a solid surface extending to the out- 
flow plane. Figure 2 shows the simplified geometry. 
At the high Mach number and high angle of attack flow 
which is being st.udiedr the pressures and densities are 
extremely lowFand the axial Mach number is predomi- 
nantly supersonicTso the influence of the simplifications 
should be localized. 
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Figure 3: Planform and frontal details of the surface 
mesh. 

The computational mesh was constructed in sev- 
eral phases. FirstT a surface mesh was generated on 
the simplified geometry using the GRIDGEN package. 21 
Based on this grid distribution and conservative esti- 
mates of the outer boundary locationr GRIDGEN was 
used to generate surface meshes on the remaining bound- 
ary faces of the computational domain (the upper and 
lower symmetry planesrthe outflow planeEand the outer 
boundary). Since all the faces of a “computational cube” 
had now been definedrthe aim briar volume grid genera- 
tor within GRIDGEN was used to generate the interior 
points. This algebriar volume grid was then smoothed 
using the Poisson solver within GRIDGENT specifying 
that grid lines should be normal to the body surface and 
symmetry planes. 

The final computational mesh is comprised of just 
over one million grid points with 121 points distributed 
stream-wise along the bodyriOl points circumferentiair 
and 91 points between the body and just outside the 
bow shock (i.e.rin the uniformrfree-stream flow). Fig- 
ure 3 shows two details (planform and frontal views) of 
the surface meshrwhile Fig. 4 depicts selected planes of 
the volume mesh (with every other grid line perpendicu- 
lar to the the body-normal direction removed for visual 
clarity). 

LAURA provides for automatic grid adaptation to 
the flowfield after an initial solution is generated. This 
automatic redistribution of points plays the delicate bal- 
ance of providing proper spacing near the wall to resolve 
the boundary layeiT clustering points in the vicinity of 
the bow shockr and also providing a minimum number 
of points outside of the bow shock structure — i. eTthat 
portion of the computational domain where freestream 
conditions prevail. This form of adaptation not only al- 
lows for a highly efficient use of the limited number of 


Figure 4: Partial view of volume grid (every other 
line perpendicular to the body-normal direction 
omitted for clarity fi 

grid points available due to computer memory restric- 
tionsrbut also provides better resolution of the strong 
bow shock by aligning the computational cells with the 
shock gradients. 


Numerical Method 

The LAURA code is a point-implicitTfinite-volume 
solver based on the upwind-biased flux difference split- 
ting of Roe. 22 The scheme utilizes Yee’s symmetric total 
variation diminishing discretization 23 to achieve second- 
order spatial accuracy while incorporating Hart.en’s en- 
tropy fix. 24 LAURA is capable of modeling both Eu- 
ler and laminar Navier- Stokes flow for a host of dif- 
ferent air chemistry assumptions: perfect gasr equi- 

libriumr chemical non-equilibriumr and thermochemical 
non-equilibrium. For this computationELAURA was pri- 
marily run in thin-layer r laminar Navier-Stokes mode 
with chemical non-equilibrium. Specific details of the 
basic algorithm can be found in Ref. 25. 

For the altitude and velocity of the flight point con- 
sidered in this st.udyr it has been shown by Hartung 
et al. R b that the flow over the Shuttle Orbiter is still 
laminar. The thin-layer assumption (neglecting viscous 
terms tangential to the body surface) is used for two 
reasons: (i) reduction of computational resources used 
(both memory and CPU time) and (ii) as outlined by 
Baldwin et a/.R' there is generally not sufficient grid res- 
olution to resolve the streamwise gradients if they were 
to become significant to the viscous energy momentum 
balance. The effect of neglecting these terms is discussed 
in a later section of the paper. 
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Temperature (K) 

Figure 5: Atomic recombination coefficients for 

Oxygen and Nitrogen as a function of temperatures, 

Aerothermodynamic Models 

Some of the major considerationsT specific to this 
flowfield application are presented below. Extensive 
details of the aerothermodynamic models used in the 
LAURA code can be found in Ref. 28. 

Finite- Rat® Chemistry 

The fiowfield chemistry is described by the kinetic 
model of Park 29 with modifications for the dissociation 
rates according to Ref. 30. HoweveiT two simplifying 
assumptions were used for this case; 

1. Collapsing the two-temperature model into a single 

temperature. 

2. Using only seven of the proposed eleven species 

The first assumption is physically valid because the 
post-shock temperatures of this flow are relatively low 
(?s6500K)rand thus the vibrationairtranslationairrota- 
tionairand electronic energy modes are near equilibrium. 
The seven species accounted for in this calculation are 
0 2 r#2r0rtVrtV0rtV0 + re-. Againrdue to the rela- 
tively low energy level of the flowrionization of nitrogen 
and oxygen would be insignificant if included. (This as- 
sumption has been re-affirmed by Hart.ung through use 
of an approximate method within Ref. 31.) 

Finite-Rate Wall Catalysis 

Zoby 32 and Scott’s 33 finite-rate wall catalysis curve 
fits were used to model the wall catalysis of the Shuttle’s 
RCG coated tiles. Figure 5 shows a plot of atomic re- 
combination coefficients with respect to temperature for 
Oxygen 32 and Nitrogen. 33 The validity ranges for both 
fits are also shown on the plot. Limited data for the 
lower temperature ranges are available from Kim 34 and 


Marinelli 35,3b which show a trend very similar to that 
given by the extrapolated curves of Zoby and Scott. The 
latter method was chosen for the present study to sim- 
plify the model and limit the computational expense of 
computing the wall catalysis. 

In addition to the above simplification of the wall 
catalysis recombination ratesT all surfaces of the shut- 
tle ( i. eTall types of reusable tile and other surface in- 
sulationrincluding the reinforced carbon-carbon (RCC) 
found on the leading edge surfaces) were assumed to fol- 
low this trend. 

Variable Wall Temperature Distribution 

The Stefan- Boltzmann relation was used to develop 
wall temperature distributions from the computed heat- 
ing rates. This relation links the emissivity of a surfacer 
er(t.he ability to emit radiation compared to a Planck 
black-body) with the heating raterf/Tof that surfacer 

q = acT 4 

where a is the Stefan-Boltzmann constant (5.67 x 
10 _8 W/m 1 2 Iv 4 ). 

The emissivities are a function of temperature and 
surface type as given by Refs. 17-19. The surface of the 
Shuttle is composed of many disparate materialsLand for 
this studyrsurfaces were considered as being either dark- 
or light-colored for determining which set of emissivity 
functions to use. 

Using these relations for emissivities and the Stefan- 
Boltzmann relationra wall • < 1 1 1 1 >< ra i nn I I I at; each cell 
face for the Shuttle surface can be determined from the 
predicted heating rater q. This iterative process is sus- 
tained until neither the predicted heating ratesr or the 
inferred wall temperature distribution are changing sig- 
nificantly (« 10° Iv). 

Flight Data 

As part of the design confirmation procedure for 
the Shuttlerthe Developmental Flight Instrumentation 
(DFI) was installed on the Orbit.er Columbia for the first 
five flights. Among these instruments were some two 
hundred thermocouples and pressure transducers dis- 
tributed over the vehicle surface. The surface temper- 
atures and heating rates for those flights are taken from 
Refs. 17-19rwhile the flight surface pressures have been 
extracted from the DFI database. 

The SILTS data represents high-spatial-resolution 
temperature measurements of the leeside surfaces. 20,3 ' 
An infrared camera was mounted in a pod at the top 
of the vertical taiirand it alternatively records the ther- 
mal environment of the aft. portion of the fuselage and 
the wing. From these images and knowledge of the Or- 
bit.er’s surface propert.iesr the surface temperatures can 
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be determined. For the SILTS flights used in this studyr 
Le.rSTS-28 and STS-32Fonly a small contingent of the 
original DFI thermocouples on the leeside of the vehicle 
remained operational. The DFI temperature and heat- 
ing data for these flights were obtained using the same 
methods as outlined in Refs. 17-19. 

Results and Discussion 

The computational resources required for this study 
are discussedrfollowed by a description of the flight con- 
ditions used for comparison. NextTa description of the 
salient features of the leeside flow is givenrfollowed by 
comparisons of surface quantities: pressuresr tempera- 
turesFand heating rates. 

Computational Resources 

Using the methods described by Weilmuenster et 
al . 38 the Shuttle Orbiter grid was cut into four stream- 
wise sections. The sections varied between 35 and 45 
cells in the streamwise directionFand a typical run would 
require approximately 70 Cray-2 hours and lOOMw of 
memory. This time includes running on the initial gridr 
aligning the gridr and converging the temperature dis- 
tribution for each block. After the four sections were 
convergedr all four were put back together into one en- 
tire grid and run for an additional 50 Cray-2 hours. For 
this latter runrthe code was altered so that it could fit 
within 125Mw of memory at the expense of efficiency. 

Flight Conditions 

A flight trajectory point was chosen so that both 
DFI data (primarily STS flights 2T 3T and 5) and the 
SILTS data (STS flights 28 and 32) would be available 
for comparison. The point is fairly early in the trajectory 
so that the flow over the Shuttle Orbiter is still laminarr 
but yet low enough (ps 70km) for the continuum assump- 
tion to remain plausible (Knudsen numbers ps 0.05). The 
freestream velocities and densities were matched for five 
flights: STS-2r STS-3r STS-5r STS-28r and STS-32 . 
Table 1 shows the freestream conditions for the matched 
trajectory point. The freestream conditions of STS-28 
were used for the numerical simulation of this study. 

Flow Structure 

Figure 6 shows an computer-simulated oil-flow rep- 
resentation of the leeside portion of the Shuttle. The 
figure clearly depicts many of the separation and re- 
attachment lines present on the vehicle. Starting from 
the leading edge of the wing and working around to the 
top portion of the vehicle: 

• a region of very low pressure/density due to rapid 
expansion of the flow around the wing tip (z.eTthe 
wing-tip vortex)r and the corresponding cross-flow 


shock is shown by the abrupt change in flow direc- 
tion just inboard of the wing tip 

• a separation line running just behind the leading 
edge of the wing strake and continuing through the 
aft wing section caused by a standing shock trig- 
gered by the rapid expansion around the strake lead- 
ing edge (z.eTthe wing strake vortex) 

• a re-attachment/shear line located diagonally from 
the beginning of the wing strake to the leading 
portion of the Orbital Maneuvering System (OMS) 
pods 

• a stagnation region just ahead of the OMS pods 

• a separation line along the upper portion of the 
body attributed to the standing cross-flow shock 
which serves to turn the flow parallel to the sym- 
metry plane 

• a smaller separation line on the aft portion of the 
OMS pods due to the cross-flow shock 

• a re-circulation region just ahead of the canopy cre- 
ated by the canopy shock 

Figure 7 shows particle traces which originate from 
just ahead of the canopy. Note that the flow from the re- 
circulation region just ahead of the canopy (Fig. 6) leaves 
the surface and goes upTover the canopy and then spirals 
back downr wetting the top section of the aft fuselage. 
The implication of this flow structure will become evi- 
dent in the discussion of the aft fuselage heating rates. 

Surface Pressures 

As shown in Fig. 8T comparisons were made with 
DFI flight data for the windward and leeward centerline 
pressures. The windward centerline pressures are in ex- 
cellent agreement. Towards the aft end of the vehicle 
(from approximately X/L=0.93Faft) the computational 
model does not accurately model the rear portion of 
the vehicleFand thus the dramatic drop and consequent 
rise in pressure is not physical. The leeward centerline 
pressures show show good agreement with the limited 
amount of flight data. 

Figure 9 shows the pressure distribution at the sixty 
percent semi-span wing station (2Y/B=0.6) as a func- 
tion of non-dimensionalized chord position (X/C). The 
elevon hinge-line occurs at X/C=0.7 for this semi-span 
stationFand the corresponding deviation can be seen in 
the predicted lower surface pressures due to the non- 
physical modeling of the aft portion of the wing. The 
leeward pressures compare favorably until approximately 
the 50 percent chord position. At this point the flight 
data shows considerably higher pressures. The effect of 
the elevon deflection might be suspectedrbut for STS-3 


5 



Table 1: Matched Flight Conditions 


Quantity 

units 

STS-2 

STS-3 

STS-5 

STS-28 

STS-32 

UTC 

sec 

75620 

56588 

51120 

47716 

33166 

h 

km 

69.2 

72.9 

72.4 

73.2 

72.8 

M® 

- 

24.3 

22.2 

23.6 

24.3 

24.4 

F® 

km/sec 

6.92 

7.01 

6.83 

7.01 

7.01 

Too 

K 

202 

248 

209 

210 

207 

P 

1 OO 

Pa 

3.33 

3.76 

3.35 

3.40 

3.32 

Poo 

kg/m 2 xlO -5 

5.75 

5.65 

5.59 

5.28 

5.59 

a 

deg 

39.4 

39.5 

40.3 

39.4 

40.0 

P 

deg 

-0.3 

-0.1 

0.2 

n/a 

n/a 

bank 

deg 

57 

64 

-64 

-66 

-71 



Figure 6: Computer simulated oil-flow on the leeside of the Shuttle Orbiter. 
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Figure 7: Streamlines originating just ahead of the canopy. 




Figure 8: Windward and leeward centerline pressure 
comparisons. 


Figure 9: Pressure distribution comparison for wing 
station 2Y/B=0.60 as a function of chord. 


and STS-5 the elevens were deflected differently (down- 
ward 1.7° and 3. 7° Respectively), This discrepancy be- 
tween predicted and measured pressures is apparently 
due to bleed-through from the gap between the inboard 
and outboard elevon. This gap is approximately 18 cm 
(7 inches) wide and occurs just inboard of this wing 
semi-span station. Since this gap was not modeled for 
the computational solutioiiPthe corresponding increase 
in pressure due to bleed-through from the windside is 
not present. 


Figure 10 shows the pressure distribution at the ten 
percent fuselage station (X/L=0.1). The meridional an- 
gleF i> I is measured from the leeward symmetry plane 
around to the windward symmetry plane. The figure 
also includes a half-section depiction of the fuselage at. 
the X/L=0.1 station. At this locationr the predicted 
and measured pressures agree very well. Again due to 
sparse pressure data for the leesideFno direct compari- 
son can be made for the top of the fuselage. The trend 
of the predicted pressures] however I is as expected for a 
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Figure 10: Cross-sectional pressure comparison at 
fuselage station X/L=0.10. 


blunt-nosed vehicle at moderate angle of attackrshowing 
a rapid expansion and subsequent re-compression while 
approaching the symmetry plane. 

Surface Temperatures 

Figure 11 is a side-by-side comparison of the lee- 
side wall temperature distribution for the aft portion of 
the Shuttle. The computed results from LAURA are 
shown on the left-hand side while the measured temper- 
atures from the STS-32 SILTS flight are shown on the 
right. The SILTS image for STS-28 flight point is nearly 
identical to STS-32 result.sr so only the one compari- 
son is presented. Recall that the LAURA temperatures 
are determined from the radiative equilibrium assump- 
tion discussed in a previous section. This equilibrium 
assumption neglects the effects of the time dependent 
terms such as conduction. 

The figure shows a band of high temperature along 
the leading edge of the wing as expected. The spotted 
nature along the forward section of the wing shown in 
the SILTS image is due to temperature variations across 
each of the hollow RCC segments which comprise the 
leading edge. This phenomena is largely due to the re- 
radiation from the windward surface of the hollow RCC 
leading edge sections through to the leeward surface (see 
Ref. 37). Also depicted by both methods is a streamwise 
streak of high temperature at approximately 60 percent 
of the wing semi-span. This corresponds to the wing-bow 
shock interaction. In additionrboth images show an ar- 
eas of low temperature located just under the wing tip 
vort.exrbelow the inboard wing-st.ra.ke vortex structurer 
and just ahead of the of the OMS pods. The SILTS im- 
age shows a streak of high heating where the windward 
flow comes through the elevon-elbvon gap at the trailing 
edge of the wing. The computational shuttle model did 


not model this gap. The rise in heating on the leading 
edge of the OMS pods is clearly depicted by both meth- 
ods. Also of interest is that even though the tail was 
neglected from the computational model of the Shuttle 
OrbiteiT the temperature distributions are remarkably 
similar on the top of the fuselagerjust. ahead of where 
the tail is located. 

Figure 12 shows the wing upper surface non- 
dimensiona.lized temperature distribution as a function 
of non-dimensiona.lized chord position for three semi- 
span stations (2Y/B): 0.40r0.60Fa.nd 0.75. Each figure 
includes DFI thermocouple dataF SILTS dataF and the 
temperature profile predicted by LAURA. The SILTS 
data for these plots was extracted from the image data 
shown in Fig. 11. The extracted SILTS data had to 
be stretched to properly register the leading and trailing 
edges. This is due to uncertainties in the optical transfor- 
mations used to obtain the planform image as discussed 
in Ref. 37. In additionr the reader is cautioned that 
the aft. portion of the LAURA solutions are not physical 
since the wing was not modeled accurately beyond the 
elevon hinge-line. 

The forty percent semi-span station (Fig. 12a ) 
shows that both sets of flight data show lower tempera- 
tures than the LAURA prediction. This might be dm In 
conduction along the RCC leading edge which was not 
computationally modeled. At the sixty percent semi- 
span station (Fig. 12b) the temperatures predicted by 
LAURA and the flight measurements are in very good 
agreement considering this semi-span station is very near 
the wing-bow shock interaction zone. The DFI ther- 
mocouple near the leading edge (X/C=0.15)rhoweveiT 
is considerably lower than either SILTS or LAURA. At 
the seventy-five percent semi-span station (Fig. 12c) the 
LAURA predictions are within five percent of both sets 
of flight data. 

Surface Heating 

Figure 13 shows heating comparisons for the wind- 
ward and leeward centerlines. The windward centerline 
values are in very good agreement with all three DFI 
flights. One isolated discrepancy occurs at the fifteen 
percent fuselage station (X/L=0.15) where the heating 
rates shown by the DFI data are non-monot.onic for a 
short stretch of the vehicle surface. This is attributed 
to an acoustic microphone which melted on one of the 
first flights and contaminated a localized downstream 
area. It should be noted that several of the tiles on all 
flights (particularly STS-5) had catalytically enhanced 
surfaces. 39 HoweveiTa.s discussed in Ref. 39rt.he area, of 
influence appears to be very localized high heating rates. 
(These offending points have been omitted for the sake of 
clarity.) The predicted leeward centerline heating rates 
show a. consistently lower trend than the measured flight. 
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Figure 11: Leeside wall temperature comparison for the aft end of the Shuttle Orbiter. 


dataP and once past X/L=0.55 not even the trend in 
heating is captured correctly. Of particular note is the 
peak heating caused by the canopy shock is predicted 
significantly lower than that measured by flight data. 
This region is just aft of a bank of four Reaction Con- 
trol System (RCS) jets which were not modeled. This 
region also contains the recessed windows and their as- 
sociated gap sealsr etc — none of which are modeled in 
detail. The low predicted heating in this area is of con- 
cern because the vortex which originates from just ahead 
of the canopy deposits its flow along the aft end of the 
fuselage (refer to Fig. 7). Since there are not very many 
grid points in the vicinity of this region (Fig. 3)rthe flow 
is not adequately resolvedrand as a consequence of this 
vortical structureFthe flow which later impinges on the 
aft end of the fuselage does not have as high of energy 
as it should. Another plausible factor contributing to 
the lower predicted heating rates is the solar radiation 
present at this altitude. Throckmorton 40 shows that so- 
lar radiation contribution to the heating rates on isolated 
leeside gauges could be up to 60 percent of the total for 
a given inclination to the Sun during reentry. 

The heating comparisons at several fuselage cross- 


section stations are shown in Figure 14. As with the 
pressuresrt.he meridional anglerTris measured from the 
leeward symmetry plane around to the windward sym- 
metry planeFand each figure includes a half-section de- 
piction of the fuselage at the particular station. 

Figure 14a shows the first fuselage station 
(X/L=0.1). This station is located midway between the 
cluster of four RCS jets and the beginning of the canopy. 
The predictedrieeward heating rates are somewhat lower 
than the measured flight data. HoweveiT as one ap- 
proaches the windward surfaces the predicted heating 
rates show excellent agreement with the flight data. 

Figure 14b represents the X/L=0.25 fuselage station 
which is located just aft of the forward facing windowsT 
but still encompasses the swept-back portion of the fuse- 
lage which accommodates the side-facing windows. This 
station shows a trend similar to the previous station: 
under-prediction when nearing the leeward centerlineF 
but excellent, agreement on the sidewall and windward 
surfaces. At this stationrhoweveiTthe trend of near the 
leeside centerline is not even captured. 

Figure 14c (X/L=0.3) is located just aft of the pre- 
vious stationTaud no longer encompasses the swept-back 
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Figure 12: Temperature distribution comparisons at 
various wing semi-span stations as a function of chord. 

portion of the fuselage due to thf side facing windows. 
Once againrthe windward heating rates agree wellTwhile 
the predicted leeside heating rates are lower than those 
measured in flight. Againrthe trend of the heating rates 



Figure 13: Windward and leeward centerline heating 
rate comparisons. 

near the leeside centerline is not properly captured. 

Figure 14d depicts the fuselage station at X/L=0.4. 
At this pointTbot.h the windward and leeward heating 
predictions agree very well with the measured flight data. 
The predicted heating near the leeside centerline is still 
slightly lower than the flight da.t.aFbut. the proper trend 
is now evident. 

The next fuselage station (X/L=0.5) is shown by 
Fig. 14e. As with the previous cross-sectional stationr 
the leeside heating rates are in slight disagreementTbut 
all the trends are captured. A slight disparity in the 
flight data is noted at a $ of approximately 75° . 

The next three sub-figures: 14fT14grand 14h show 
the Computed and measured heating rates for fuselage 
stations of X/L=0.6r 0.7r and 0.8. All show similar 
trends: excellent agreement on the windward and side 
surfacesrbut. over-prediction of the heating rates on the 
upper surface of the wing and under-prediction on the 
lee surface of the fuselage. 

General Observations 

As a check of the thin-layer Navier-St.okes a.ssump- 
tionrthe neglected viscous terms are added in two stages. 
The first, stage is the addition of the circumferential com- 
ponent.. The inclusion of this term involves no additional 
memory and is relatively inexpensive to compute. With 
this additional t.errnr the changes in the heating rates 
are highly localized and generally less than five percent.. 
These occur in the wing/bow-shock interaction regionr 
in the vicinity of the wing-tip cross-flow shockra.nd just, 
ahead of the OMS pods. The second stage is to use the 
full Navier-St.okes equations. This increases the code’s 
memory requirement, by approximately fifteen percentT 
but. the computational overhead is no more than five 
percent.. Although the arguments of Baldwin et a/.P' 
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Figure 14: Cross-sectional heating comparisons at various fuselage stations along the vehicle. 
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Figure 14: 

indicate thatT in generair the effects of t.hgb additional 
terms would be negligibler results demonstrate that in 
confined regionsrt.he leeside heating rates change by as 
much as thirty percent. For example the heating rates 
on the sidewall of the fuselage in the region below the 
shear-line extending between the start of wing-st.rake (re- 
fer to Fig. 6) decrease significantly. This is as might 
be expected since this is a highly vortical region of the 
flow. The heating increases by twenty percent along both 
cross-flow shocks on the wing. AlsoFa decrease in heat- 
ing is noted at the front of the OMS pods and along the 
lower OMS pods/fuselage juncture. 

The procedures used to calculate the heating rates 
from the DFI data 1 1-19 make assumptions concerning 
the integrity of the thermocouple/thermal protection 
material joint. For the flexible materials used through- 
out the leeside of the Shuttle Orbit.eiT this is not nec- 
essarily a valid assumption. Estimates of error in the 
procedures outlined in Refs. 17-19 do not include the 
errors present in the thermocouple measurements them- 
selves. Thereforer considerable unknowns are present 
in the thermocouple flight data and corresponding heat 
transfer rates. 

The atomic recombination coefficients (i.e.r wall 
catalysis) for the surfaces of the Shuttle Orbit.er are very 
difficult to measure precisely. This is evident in the wide 
spread of atomic rate coefficients determined experimen- 
tally for only one type of surface coating material in a 
limited temperature range. 41 As mentioned earlieiT all 
the surfaces of the Shuttle Orbit.er were modeled as hav- 
ing the same catalycity function with temperature. This 
is definitely not the case because there are many differ- 
ent types of surface coatings. For examplert.he Advanced 
Flexible Reusable Surface Insulation (AFRSI) which is 



Concluded. 

used in many regions on i In- the leeside of the Orbit.er is 
not coated with the same RCG coating that, the catalyc- 
ity model considers. This global model is also not cor- 
rect. for the RCC found on the wing leading edges and 
the fuselage nose cone. Other factors such as surface ag- 
ing and contamination are completely unknownrbut. are 
expected to produce only second-order effects. 

Conclusions 

A computational solution was presented for a. reen- 
try condition of the Shuttle Orbit.er at. 40° angle of attack 
for an altitude of 70 km. The initial results demonstrate 
that, the flow structure is quite complexr maintaining 
several vortical struct.uresrsepara.t.ion linesFshock-shock 
intera.ct.ionra.nd regions of low densit.yrfrozen flow. 

The computed pressures agree well with flight data, 
on the windward and leeward centerlines and for a. fuse- 
lage cross-section. The pressures along a. wing serni- 
spa.n station also show good agreement, for the wind- 
siderbut. differ considerably on the leeside near the in- 
board/outboard elevon gap which was not modeled com- 
putationally. 

The temperature distributions on the leeside of the 
wing are well predicted by LAURA. The wing-bow shock 
interaction region is evident as well as the peak hatting 
on the wing leading edge. This offers confidence in ac- 
curately predicting the leeside thermal environment for 
design considerations. 

Using the thin-layer Na.vier-St.okes assumption for 
the leeside flow appears to be adequate for the flight con- 
dition studied in this paper. The effects of adding the 
remaining viscous terms produced only localized changes 
in the heating rates for the leeside flowr and had abso- 
lutely no effect, on the windside of the vehicle. Of course 
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a grid refinement study is necessary to completely an- 
swer this question. 

The heating ratesr overalir compare very favorably 
with those measured by DFI flight data. The heating 
rates predicted by LAURA for the top of the aft fuse- 
lage are considerably lower than that measured by DFI 
data. It is expected that better resolution of a vorti- 
cal structure emanating from the canopy region might 
rectify this. 
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